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ABSTRACT 

Experimental  investigations  have  been  conducted  to 
examine  the  general  problem  of  convective  heat  transfer  in  a 
rectangular  cross  section  channel  with  both  straight  and 
curved  passages.   The  experimental  configuration  was  modeled 
as  infinite  parallel  plates  with  one  wall  at  a  constant  heat 
flux  and  the  opposite  wall  adiabatic.   The  air  flow  through 
the  channel  was  heated  through  the  Joulean  heating  of  a  wall 
composed  of  Temsheet.   Liquid  crystals  applied  to  the  surface 
of  the  Temsheet  were  used  to  confirm  the  presence  of  the 
Taylor-Gortler  vortices .   Average  Nusselt  number  measurements 
were  used  to  investigate  their  effect  on  heat  transfer  from 
the  curved  test  section.   Experiments  were  conducted  for 
steady  state  laminar  flow  that  was  hydrodynamically  developed 
and  thermally  developing.   It  was  confirmed  that  the  Taylor- 
Gortler  vortices  enhanced  the  heat  transfer  process  in  the 
curved  section.   The  experimental  data  was  compared  with 
several  analytic  studies  with  good  results . 
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I.   INTRODUCTION 

A.   DESCRIPTION  OF  TAYLOR-GORTLER  VORTICES 

A  considerable  amount  of  data  and  indirect  evidence,  [l, 
2,  3]  among  others,  has  shown  that  the  fully  developed  laminar 
flow  of  viscous  fluids  along  a  concave  wall  does  not  remain 
two  dimensional.   Instead,  the  flow  forms  spiral  vortices  in 
counter-rotating  pairs  which  frequently  possess  a  regularly 
spaced  cellular  structure.   This  phenomenon,  known  as 
Taylor-Gortler  vortices,  occurs  as  a  result  of  flow  insta- 
bility induced  by  variations  in  the  centrifugal  forces  acting 
on  fluid  particles.   The  Taylor-Gortler  vortices  are  laminar 
in  nature  with  their  axes  aligned  in  the  direction  of  the 
principal  fluid  flow  and  with  secondary  flow  velocities  in 
both  perpendicular  directions .   An  example  of  this  flow 
pattern  is  shown  in  Figure  1. 

In  a  channel  that  is  curved  in  the  streamwise  direction, 
such  as  that  used  in  this  study,  the  fluid  particles  near 
the  center  are  subjected  to  higher  centrifugal  forces  than 
those  experienced  by  the  slower  moving  fluid  particles  near 
the  boundary  of  the  concave  wall.   Consequently,  the  tendency 
is  for  the  fluid  particles  in  the  center  of  the  channel  to 
move  outward  to  the  concave  wall.   Under  this  influence, 
particles  near  the  wall  are  displaced  initially  in  the 
spanwise  direction  and  subsequently  continue  radially  inward, 

12 


Figure  1.  Taylor-Gortler  Vortex  Flow  Pattern 


13 


replacing  the  central  fluid  particles.   Once  in  the  center 
of  the  channel  the  fluid  particle  will  come  under  the  in- 
fluence of  the  higher  free  stream  velocity  and  resultant 
centrifugal  forces.   This  action  forces  movement  back  toward 
the  concave  wall  causing  the  entire  cyclic  rotation  process 
to  propagate . 

Although  the  Taylor-Gortler  vortices  are  laminar,  it  is 
quite  possible  that  they  play  an  important  role  in  the  transi- 
tion from  the  laminar  to  the  turbulent  flow  regime.  [2] 
Additionally,  the  presence  of  the  vortices  have  been  thought 
to  account  for  the  increase  in  heat  transfer  from  curved 
walls  due  to  mass  transfer.  [4] 

The  Taylor-Gortler  vortices  also  have  similarities  to 
other  vortex  flow  patterns  such  as  the  longitudinal  vortex 
rolls  developed  in  the  laminar  forced  convection  heating 
of  fluid  layers  between  parallel  plates.  [5]   In  another 
field  it  has  further  been  hypothesized  that  the  cross-hatching 
phenomenon  observed  in  re-entry  vehicles  can  be  explained 
by  the  presence  of  centrif ugally  induced  secondary  motion 
characteristic  of  the  Taylor-Gortler  vortices.  [6] 

It  is  apparent  that  the  understanding  of  these  vortices 
and  their  effect  on  heat  transfer  and  fluid  flow  character- 
istics could  significantly  effect  many  areas  of  fluid 
mechanics.   Recently  it  has  been  proposed  [7,  3,  9  J  that 
concave  wall  curvature  can  have  an  advantageous  effect  on 
film  cooling  over  turbine  blades. 
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Many  applications  of  fluid  mechanics  and  heat  transfer 
involve  curved  surfaces  or  curved  channels  and  the  under- 
standing of  the  fluid  mechanics  in  this  geometry  will  greatly 
enhance  engineering  designs  and  capabilities. 

B.   BRIEF  HISTORY 

Instability  of  an  inviscid  fluid  in  a  curved  flow  was 
first  considered  by  Lord  Raleigh  in  1916.  [10  ]   He  concluded 
that  the  stability  criterion  for  an  inviscid  fluid  had  as  a 
necessary  condition  that  the  circulation  always  increased 
with  the  radius.   In  essence,  this  meant  that  the  product 
of  the  local  circumferential  velocity  and  the  corresponding 
local  radius  of  curvature  cannot  decrease  with  increasing 
radial  distance . 

In  1923,  G.  I.  Taylor  expanded  the  scope  of  Raleigh's 
stability  criterion  to  include  viscous  fluids,  [l ,  11]   His 
extensive  analytic  and  experimental  investigation  focused  on 
the  flow  of  viscous  fluids  between  a  rotating  inner  cylin- 
der and  a  stationary  outer  cylinder.   Taylor  ascertained  that 
such  couette  flows  become  unstable  when  a  characteristic, 
non-dimensional  parameter  exceeded  a  critical  value  of  M-l. 
Referred  to  as  the  Taylor  number,  this  parameter  was  defined 
as: 


ia    v  /  R.      Ke   R. 
l  l 
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In  this  relationship  d  represents  the  spacing  between  the 
concentric  cylinders,  assumed  small  when  compared  to  the 
radius  of  curvature,  R..   Re  is  the  Reynolds  number  based  on 
U,  the  circumferential  velocity  of  the  inner  cylinder  and  d, 
the  annulus  width .   When  flow  conditions  were  such  that  the 
Taylor  number  exceeded  this  critical  value,  the  secondary 
flow  observed  in  the  gap  between  the  stationary  outer  cylin- 
der and  the  concentrically  rotating  inner  cylinder  was  in  the 
form  of  longitudinal  vortices  whose  axes  were  aligned  in  the 
circumferential  direction.   These  vortices  were  given  the  name 
Taylor  vortices. 

A  similar  form  of  instability  occurs  when  a  viscous 
fluid  flows  in  a  curved  channel  due  to  a  pressure  gradient 
acting  along  the  concave  wall.   W.  R.  Dean  f 12 J  first  addressed 
this  problem  analytically  in  1928.   Specifically,  he  consi- 
dered flow  in  a  channel  formed  by  two  concentric  cylinders  with 
the  restriction  that  the  cylinder  gap  distance  be  small  com- 
pared to  the  radius  of  the  inner  cylinder.   Dean  concluded 
that  the  flow  instability  would  initiate,  and  the  secondary 
flow  vortices  (similar  to  the  Taylor  vortices)  would  form, 
when  the  Dean  number  was  greater  than  36.   The  Dean  number 
was  then  defined  as: 

De  =  Re  /  ^r- 

Ri 
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In  this  case,  the  Reynolds  number  was  based  on  the  mean  velo- 
city of  the  undisturbed  flow  and  d  represented  the  channel 
half  width.   The  analytical  work  by  Dean  was  later  verified 
by  W.  H.  Reid  [13]  using  an  approximate  numerical  solution. 

In  1940,  H.  Gortler  [2]  studied  the  influence  of  small 
disturbances  on  the  stability  of  laminar  boundary  layer 
profiles.   Gortler  found  that  these  disturbances  closely 
resembled  those  described  by  G .  I.  Taylor  which  led  to  fluid 
flow  instability  in  the  form  of  vortices.   In  his  numerical 
calculations,  Gortler  concluded  that  the  amplified  distur- 
bances were  produced  only  on  the  concave  walls;  nonetheless, 
the  overall  flow  regime  seemed  to  remain  laminar.   The  appro- 
ximate analytic  results  that  Gbrtler  had  developed  were 
initially  verified  with  an  exact  solution  by  G.  Hammerlin  as 
reported  by  H.  Schlicting  [14],  and  later  substantiated  by 
an  extensive  numerical  analysis  by  A.  M.  0.  Smith  [3].   M. 
Kelleher,  et  al.,  [15]  and  S.  Winoto ,  et  al.,  [16],  among 
others,  have  been  able  to  verify  these  numerical  predictions 
by  using  hot  wire  anemometry,  laser  Doppler  systems  and  flow 
visualization  techniques  to  describe  the  motion  of  the  secon- 
dary flow.   A  non-linear  analysis  of  Gortler  vortices  was 
conducted  by  Y.  Aihara  [17]  in  1976. 

Growing  interest  in  the  various  effects  of  the  secondary 
flow  associated  with  Taylor-Gortler  vortices  led  to  more 
specific  studies  concerning  the  vortice  influence  on  heat 
transfer  in  both  turbulent  and  laminar  flow.   In  1955, 
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F.  Kreith  [4]  reported  that  the  heat  transfer  from  a  heated 
concave  wall  was  considerably  higher  than  that  from  a  convex 
wall  of  identical  curvature  under  similar  conditions  of  tur- 
bulent flow. 

In  the  special  cases  of  very  high  or  very  low  Prandtl 
number,  L.  Persen  [18]  related  the  increased  heat  transfer 
from  curved  walls  to  the  presence  of  Taylor-Gortler  vortices. 
He  determined  that  their  overall  effect  was  to  enhance  heat 
transfer  through  the  boundary  layer. 

In  1970,  McCormack,  et  al . ,  [19]  published  the  results 
of  their  initial  experimental  work  regarding  the  effect  of 
secondary  flow  on  heat  transfer.   From  the  data  obtained  they 
concluded  that  a  complete  theoretical  explanation  of  the  heat 
transfer  effects  required  that  all  non-linear  terms  be  re- 
tained in  the  flow  equations.   R.  Kahawita  and  R.  Meroney  [20] 
considered  this  hypothesis  and  further  reported  that  both 
the  high  order  terms  and  the  normal  velocity  components  of 
the  main  flow  were  essential  to  any  calculations  at  small 
wave  numbers . 

In  1966,  Y.  Mori  and  Y.  Uchida  [5]  considered  the  fully 
developed  forced  convection  heat  transfer  between  horizontal 
plates.   In  the  results  of  their  investigation  they  reported 
that  as  the  temperature  difference  between  the  plates  was 
increased  a  critical  value  developed  beyond  which  longitudinal 
vortex  rolls,  similar  to  the  Taylor-Gortler  vortices,  would 
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form  with  their  axes  in  the  streamwise  direction.   M.  Akiyama, 
et  al.,  [21]  confirmed  the  presence  of  these  longitudinal 
vortices  at  a  critical  Rayleigh  number  of  1708  in  cases  where 
the  plate  is  heated  from  below. 

Although  heat  transfer  in  curved  duct  flow  has,  and  will 
continue  to  be  investigated  experimentally,  considerable 
emphasis  has  been  directed  recently  towards  heat  transfer 
predictions  using  numerical  approximation  for  flows  in  curved 
duct  geometries.   In  1970,  K.  Cheng  and  M.  Akiyama  [22] 
developed  a  numerical  solution  for  forced  convection  heat 
transfer  in  laminar  curved  rectangular  channels  by  using  a 
point  successive  over-relaxation  method.   For  similar  flow 
and  boundary  conditions,  G.  Yee  and  J.  A.  C.  Humphrey  [2  3] 
used  the  fully  elliptic  forms  of  the  transport  equations  to 
compose  a  numerical  solution.   However,  both  of  these  numer- 
ical solutions  were  limited  to  small  aspect  ratios  of  essen- 
tially square  cross  section. 

In  contrast  to  the  large  volume  of  published  numerical 
approximations  to  this  heat  transfer  problem,  the  number  of 
experimental  and  analytical  investigations  have  been  very 
limited.   Results  for  a  fully  developed,  constant  wall  heat 
flux,  curved,  square  cross  section  channel  were  reported  by 
Y.  Mori,  et  al .  [24],  and  results  for  larger  aspect  ratios 
were  reported  by  M.  Durao .  [25] 
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II.   INTENT  OF  THE  STUDY 

In  view  of  the  possible  engineering  design  implications 
of  the  Taylor-Gortler  vortex  phenomenon,  the  primary  purpose 
of  this  study  was  to  investigate  the  effects  of  the  secondary 
flow  on  the  heat  transfer  rates  in  a  curved  channel  of  rec- 
tangular cross  section. 

It  was  expected  that  the  presence  of  the  Taylor-Gortler 
vortices  would  result  in  the  secondary  flow  velocity  com- 
ponents transporting  the  heated  fluid  particle  from  the  cur- 
ved wall  to  the  center  of  the  flow  channel.   This  fluid 
motion  would  displace  the  cooler  fluid  particle  in  the  center 
of  the  channel  causing  it  to  move  toward  the  heated  curved 
wall.   It  was  anticipated  that  this  cyclic  motion  would 
measurably  improve  the  heat  transfer  process  for  all  fluid 
velocities  in  the  laminar  flow  regime.   Secondary  investiga- 
tions were  conducted  to  determine  the  additional  influence 
of  buoyant  effects  on  the  development  of  heat  transfer  in 
the  curved  channel.   It  was  expected  that  the  heat  transfer 
rates  would  be  greater  when  the  buoyant  forces  were  in  the 
direction  of  the  main  flow. 

The  investigation  was  conducted  using  a  rectangular  cross 
section  channel  that  contained  both  a  straight  and  a  curved 
test  section.   The  results  obtained  in  the  straight 
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section  served  as  a  baseline  and  as  such  were  compared  to 
the  data  obtained  in  the  curved  section.   The  complete 
process  was  accomplished  for  two  distinct  channel  orienta- 
tions:  (1)  buoyant  forces  aligned  with,  and  (2)  buoyant 
forces  opposed  to,  the  main  direction  of  the  fluid  flow. 

Straight  section  data  were  compared  with  analytic  solu- 
tions provided  by  P.  McCuen,  et  al . ,  [26],  P.  M.  Wors^e- 
Schmidt  [27]  and  R.  K.  Shah  and  A.  L.  London.  [28]   The 
buoyant  effects  were  qualitatively  compared  with  the  data 
contained  in  a  numerical  computation  by  R.  Chilukuri  and 
J.  A.  C.  Humphrey.  [29] 
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III.   EXPERIMENTAL  WORK 

A.   DESCRIPTION  OF  THE  APPARATUS 

A  rectangular  cross  section  channel  built  as  described 
in  reference  25  was  used  to  achieve  the  objectives  of  this 
study.   The  channel  was  manufactured  from  two  0.6  35  cm  thick 
sheets  of  Plexiglas  separated  by  0.6  35  cm  thick  spacers  whose 
inside  lengths  served  as  the  sides  of  the  channel. 

As  shown  in  Figure  2,  the  channel  was  composed  of  a 
straight  portion  12  2.0  cm  long,  an  intermediate  curved 
section  which  formed  a  180  degree  arc  and  a  final  short 
straight  section.   The  radius  of  curvature  of  the  interior 
concave  wall  of  the  curved  portion  was  30.5  cm.   The  channel 
was  0.6  35  cm  high  and  25.  M-  cm  wide  with  an  aspect  ratio  of 
40  and  a  cross  sectional  area  of  16.13  square  centimeters. 
A  cross  sectional  view  of  the  channel  is  shown  in  Figure  3 . 

The  working  fluid  was  air  at  room  temperature  which 
entered  the  flow  channel  through  an  entrance  bell  constructed 
of  Plexiglas  .    The  entrance  bell  was  designed  and  manufactured 
in  accordance  with  ASME  nozzle  standards  with  an  elliptic 
curve  based  on  a  ten  inch  major  and  one  inch  minor  axis.   A 
permeable  cloth  covered  the  entrance  nozzle  to  reduce  the 
introduction  of  dust  and  other  foreign  matter  into  the  flow 
channel.   The  flow  at  the  end  of  the  channel  exited  through 
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an  aluminum  exhaust  bell  connected  by  flexible  tubing  to  a 
Fisher  and  Porter  Company  variable  area  flow  meter,  model 
number  10A35  65A.   The  rotameter  had  a  10  0%  full  scale  rating 
of  0.314  cubic  meters  of  air  per  minute  (11.1  standard  cubic 
feet  of  air  per  minute)  at  one  atmosphere  and  70°  F. 

The  flow  of  air  was  drawn  through  the  channel,  the  flexi- 
ble tubing  and  the  rotameter  by  a  single  electrically  driven, 
Cadillac  centrifugal  blower,  model  6-12.   The  blower  speed 
was  controlled  by  varying  the  motor  voltage  with  a  General 
Radio  Company  Variac  Autotransformer  type  W10MT3. 

For  the  purpose  of  obtaining  experimental  heat  transfer 
data,  two  test  sections  were  constructed;  (1)  a  straight 
test  section,  29.2  cm  in  length,  located  in  the  straight 
portion  of  the  channel  downstream  of  the  hydrodynamic  entrance 
region,  and  (2)  a  curved  test  section  28.3  cm  in  length,  sub- 
tending a  circular  arc  of  5  3.1°  in  the  lower  quadrant  of  the 
curved  portion  of  the  channel.   The  area  of  the  straight 
test  section  was  741.9  square  centimeters  and  that  of  the 
curved  test  section  was  717.7  square  centimeters. 

In  each  test  section  the  outer  wall  was  detached  and 
modified  as  shown  in  Figures  H    and  5.   Temsheet,  which  is  a 
carbon  impregnated  porous  paper  with  the  property  of  uniform 
electrical  resistivity,  was  glued  to  the  entire  interior 
surface  of  the  outer  wall  in  each  section.   The  flow  of  air 
could  then  be  heated  through  the  Joulean  heating  of  the  Tem- 
sheet.  Since  the  Temsheet  had  an  electrical  resistance  that 
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was  slightly  temperature  dependent,  a  precision  resistor  with 
an  electrical  resistance  of  2.0236  ohms  was  connected  in 
series  with  the  Temsheet  to  facilitate  calculating  the 
instantaneous  power  supplied.   A  LAMBDA  Regulated  D.C.  Power 
Supply  model  LK3M-5A  FM  supplied  electrical  current  to  heat 
the  Temsheet  to  desired  levels. 

The  variables  to  be  measured  in  the  investigation  were: 
the  flow  inlet  temperature  at  the  entrance  to  each  test 
section  CTT„) ,  the  flow  outlet  temperature  at  the  exit  of 
each  test  section  (TnnT) ,  the  wall  temperature  of  the  heated 
plate  in  each  test  section  (T   ) ,  the  wall  temperature  oppo- 
site the  heated  plate  in  each  test  section  (T  . ) ,  the  flow 

wi 

exit  temperature  at  the  end  of  the  channel  (T„yTT) ,  the  temper- 
atures between  the  two  insulation  layers  for  each  test  section 
(TT,TC),  and  the  voltages  across  the  precision  resistor  (V-.^) 

and  the  heater  (Vu).   Copper-Constantan  thermocouples  were 

rl 

used  to  measure  all  the  above  temperatures. 

Liquid  Crystals 

Cholesteric  Liquid  Crystals  were  employed  to  provide 
qualitative  measurements  of  the  surface  temperature  distri- 
bution on  the  Temsheet .   The  liquid  crystals  exhibit  drama- 
tic changes  in  color  for  a  small  temperature  differential. 
The  liquid  crystals  would  allow  a  reasonable  qualitative 
comparison  to  be  made  between  the  heat  transfer  processes  on 
each  test  section  by  observation  of  the  color  pattern 
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displayed  for  various  temperature  gradients  and  flow  velo- 
cities.  Specifically,  the  occurrence  of  multi-color  stripes 
aligned  in  the  flow  direction  on  the  curved  wall  would  indi- 
cate regions  of  high  and  low  temperatures  and  confirm  the 
presence  of  Taylor-Gortler  vortices. 

To  select  the  appropriate  liquid  crystal  system  for 
application  to  the  surface  of  the  Temsheet,  consistent 
values  of  flow  inlet  temperature,  Reynolds  number,  Prandtl 
number  and  Nusselt  number  were  assumed  and  an  approximate 
energy  balance  was  developed: 

Q  =  4  Cpair(T0UT  "  TIN>  =  E  APL  AT 


or 


AT      .  Re  Pr    (Ac} 


T0UT  "  TIN     Nu     (APL) 


where 


T    +  T 
AT  =  T    -   IN    x0UT 
wo        2 


A  desired  ten  degree  Celsius  difference  between  the 
inlet  and  outlet  temperatures  required  an  average  wall 
temperature  of  forty-five  degrees  Celsius.   Taking  into 
consideration  the  possible  losses  in  the  system,  the  required 
power  to  achieve  that  average  wall  temperature  was  estimated 
to  be  a  maximum  of  40  Watts. 
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Upon  completion  of  a  series  of  preliminary  experiments  it 
was  established  that  a  combination  of  liquid  crystals  would 
be  required  to  adequately  delineate  the  average  wall  tempera- 
ture and  the  vortex  phenomenon.   The  test  sections  were 
therefore  sprayed  with  a  mixture  of  three  NCR  Microencapsu- 
lated Liquid  Crystals  (R-41,  R-45  and  R-49)  prepared  in 
accordance  with  reference  30.   Using  an  artist's  air  brush, 
approximately  twenty  coats  of  the  liquid  crystal  mixture  were 
applied  directly  to  the  Temsheet  in  order  to  ensure  a  sharp 
color  contrast  and  definition  within  the  response  curve  of 
each  liquid  crystal.   It  was  hoped  that  the  resultant  color 
patterns  could  be  recorded  photographically. 

A  total  of  sixty-four  copper-constantan,  glass  insulated, 
30  gauge  thermocouples  were  constructed  and  placed  in  specific 
locations  to  measure  the  required  temperatures.   All  ther- 
mocouples were  permanently  wired  to  a  Hewlett-Packard  Data 
Acquisition  System,  Model  2010C,  which  provided  an  automatic 
printed  record  of  the  instantaneous  data  for  each  thermo- 
couple position.   The  full  set  of  thermocouples  were  connected 
in  series  to  a  Kaye  Instruments  Ice  Point  Reference,  Model 
K-1M-0-4.   In  each  test  section  five  thermocouples  were  inser- 
ted between  the  first  and  second  layers  of  the  insulation 
with  an  additional  five  inserted  between  the  second  and  third 
layers.   Each  set  of  five  thermocouples  were  connected  in 
parallel  to  read  an  averaged  temperature  value  (TTM„)  between 
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the  respective  insulation  layers.   A  set  of  four  thermo- 
couples, also  connected  in  parallel,  were  inserted  into  the 
flow  channel  at  each  of  the  following  locations :   the  en- 
trance to  the  channel,  the  exit  of  the  straight  test  section, 
and  the  exit  of  the  curved  test  section.   At  each  location 
the  thermocouples  were  distributed  in  a  spanwise  interval  of 
5.0  cm  and  at  a  depth  of  0.3  cm  to  allow  the  thermocouples 
to  read  a  bulk  temperature  for  the  air  flow.   In  each  test 
section  eight  thermocouples  were  placed  in  direct  contact 
with  the  Temsheet  through  very  small  diameter  holes  drilled 
in  the  Plexiglas   above  the  Temsheet  and  the  thermocouple 
beads  were  electrically  insulated  from  the  Temsheet  with 
ENMAR  Heat  Resisting  Glyceryl  Phthalate.   An  additional 
eight  thermocouple  beads  were  embedded  in  the  Plexiglas 
channel  opposite  those  positioned  in  the  Temsheet.   Figure 
6  provides  a  detailed  sketch  of  this  thermocouple  arrange- 
ment.  The  thermocouples  in  contact  with  the  Temsheet  pro- 
vided an  average  outer  wall  temDerature  (T   )  and  confirmed 
°  wo 

temperature  distributions  indicated  by  the  liquid  crystal 

patterns.   Those  in  the  Plexiglas   channel  wall  provided 

an  average  inner  wall  temperature  (T  . )  and  assisted  in 

wi 

estimating  the  losses  in  the  system  due  to  radiation.   All 
the  thermocouples  were  calibrated  using  a  ROSEMOUNT  Communta- 
ting  Bridge,  model  920A,  and  a  ROSEMOUNT  Constant  Temperature 
Bath,  Model  913A. 
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The  thermal  insulation  consisted  of  Armstrong  ARMAFLEX 
22  Sheet  Insulation  which  is  a  flexible  foam  plastic  material. 
Three  layers  of  insulation  were  used,  each  one  with  a 
thickness  of  0.635  cm  and  covering  an  area  Slightly  larger 
than  each  test  section.   The  insulation  was  fixed  in  place 
with  the  use  of  an  Adhesive  Heat  Resistant  Ventilation  Duct 
Tape  as  shown  in  Figure  7.   The  entire  length  of  the  channel 
was  sealed  with  GENERAL  ELECTRIC  Silicone  Rubber  sealant 
caulk . 

Two  aluminum  electrodes,  0.318  cm  thick,  were  inserted 
between  the  Temsheet  and  the  Plexiglas   as  shown  in  Figure  3 . 
Each  plate  consisting  of  the  Temsheet,  the  liquid  crystals, 
the  thermocouples  and  the  aluminum  electrodes  was  affixed  in 
the  channel  using  Teflon  screws  to  provide  electrical  insula- 
tion . 

B.   EXPERIMENTAL  PROCEDURES 

Experiments  were  conducted  for  each  test  section  following 
the  same  set  of  procedures.   The  data  obtained  in  the  straight 
test  section  provided  baseline  information  for  later  curved 
section  runs,  either  aligned  with  or  opposed  to  the  buoyancy 
forces .   The  volumetric  flow  rate  was  varied  incrementally 
from  0.063  to  0.220  cubic  meters  of  air  per  minute  to  corres- 
pond to  values  of  Reynolds  numbers  in  the  laminar  flow  regime. 
Data  was  recorded  and  calculations  were  performed  for  each 
flow  rate  increment. 
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The  appropriate  electrical  power  required  to  heat  the 
Temsheet  to  a  temperature  level  that  would  initiate  the 
liquid  crystal  response  was  determined  through  preliminary 
experiments . 

Estimates  of  the  time  required  for  the  experiment  to 
reach  steady  state,  decisions  on  what  data  would  be  recorded, 
and  frequencies  at  which  that  data  would  be  obtained  were  also 
based  on  preliminary  runs. 

The  determination  as  to  what  constituted  a  steady-state 
condition  required  careful  consideration  in  light  of  the  in- 
fluence that  the  multiple  input  variables  could  have  on  the 
results.   The  final  criteria  for  steady-state  was  based  on 
three  input  variables;  the  flow  rate  percentage,  the  outside 
wall  temperature  (Temsheet),  and  the  heater  voltage.   When 
these  input  conditions  had  a  variation  of  less  than  two 
percent  over  a  ten  minute  interval,  the  experiment  was  consi- 
dered to  have  reached  steady-state.   Under  these  guidelines 
the  assumptions  of  steady  flow  and  constant  heat  flux  from 
the  upper  wall  were  valid. 

The  time  required  to  reach  steady-state  was  found  to  be 
between  two  and  two  and  one-half  hours.   Accordingly,  all 
runs  were  two  and  one-half  hours  in  duration.   Data  was 
recorded  at  one  hour  intervals  for  the  first  two  hours  and  at 
ten  minute  intervals  thereafter. 
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The  instantaneous  values  for  all  the  temperatures  were 
automatically  printed  by  the  digital  recorder  in  milivolts . 
A  standard  thermocouple  table  was  used  in  conjunction  with 
previously  obtained  thermocouple  calibration  data  to  convert 
the  milivolts  to  degrees  Celsius.   The  volumetric  flow  rate 
was  read  directly  from  the  rotameter.   The  precision  resis- 
tor voltage  (Vp„)  and  the  heater  voltage  (V„)  were  recorded 

in  volts  and  the  instantaneous  voltage  supplied  (0  )  was 

-p 

calculated  using  the  following  relationship: 
V    V 

Q  =  pg  H 

P      RPR 
where  R_R  was  2.0  236  ohms,  the  measured  electrical  resistance 
of  the  precision  resistor  in  the  circuit. 

The  steady  state  response  of  the  liquid  crystals  was 
recorded  photographically  using  a  NIKOMAT  (FTN)  3  5mm  single 
reflex  camera  with  a  85  -  205  mm  zoom  lens.   The  film  was 
25  ASA  Kodachrome  with  the  shutter  speed  and  lighting  con- 
trolled manually. 

All  data  reduction  was  accomplished  by  desk-top  computer 
programs  and  stored  routines. 
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IV.   PRESENTATION  OF  DATA 

A.   ANALYSIS 

The  flow  examined  in  this  study  was  laminar  and  steady, 
thus  in  a  general  sense,  plane  Poiseuille  flow.   The  Temsheet 
which  possessed  the  property  of  uniform  electrical  resistivity 
approximated  a  constant  axial  and  peripheral  wall  heat  flux. 
Temperature  measurements  on  the  Plexiglas   wall  opposite  the 
heated  surface  indicated  negligible  heat  loss  through  that 
wall,  i.e.,  adiabatic.   With  an  aspect  ratio  of  40  the  chan- 
nel was  approximated  by  infinite  parallel  plates.   The 
straight  portion  of  the  flow  channel,  upstream  from  the 
straight  test  section,  was  of  sufficient  length  to  ensure 
that  the  flow  was  hydrodynamically  developed  for  all  flow 
velocities  in  the  laminar  range.   For  the  purposes  of  this 
study,  the  flow  was  considered  laminar  when  the  value  of  the 
Reynolds  number  was  less  than  2000.   Based  on  the  above 
assumptions,  the  experimental  configuration  was  then  modeled 
as  laminar  forced  convection  between  parallel  plates  with 
hydrodynamically  developed  and  thermally  developing  flow. 
The  boundary  conditions  are  fundamental  boundary  conditions 
of  the  second  kind  with  one  wall  at  a  specified  axial  and 
peripheral  heat  flux  and  the  other  wall  insulated.  [28] 
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The  appropriate  variables,  corresponding  equations  and  solu- 
tion procedures  are  contained  in  the  following  paragraphs. 
The  average  heat  transfer  coefficient  from  the  heated 
wall  to  the  flow  of  air  was  defined  by  a  form  of  Newton's 
law  of  cooling: 

Q  .   =  h  ADT  AT 
xair      PL 

In  this  expression,  Q  .   was  the  heat  convected  to  the  air, 
-  xair 

h  was  the  average  heat  transfer  coefficient,  Ap,  was  the 
surface  area  of  the  heated  Temsheet,  and  AT  was  the  tempera- 
ture difference  between  the  average  wall  temperature  (T   ) 

WO 

and  the  bulk  temperature  of  the  fluid  CT„) . 

The  actual  heat  convected  to  the  air  (Q  .  )  was  then 
independently  calculated  with  the  equation: 

Q  .   =  mC   .  (TnTTrp-TTV) 
xair      pair   OUT   IN 

where  C   .   was  the  specific  heat  of  air  at  constant  pressure 
pair  r  r 

and  m  was  the  mass  flow  rate  of  the  air.   The  mass  flow  rate 
was  the  product  of  the  volumetric  flow  rate  (q)  and  the  local 
flow  density  in  the  rotameter  (p).   The  volumetric  flow  rate 
was  corrected  to  standard  air  conditions  for  the  rotameter 
and  the  local  flow  density  was  calculated  assuming  perfect 
gas  behavior  at  constant  pressure  and  using  the  flow  exit 
temperature  (Tp„-j.„)  . 
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The  average  Nusselt  number  was  determined  from  the  aver- 
age heat  transfer  coefficient  using  the  following  expression 

77-    hd 
Nu  = 


K  - 

air 


where  d  represented  the  height  of  the  channel  and  K  .   was 

43  air 

the  thermal  conductivity  of  the  air  evaluated  at  the  bulk 
temperature  (Tg) . 

The  Reynolds  number  was  calculated  for  all  the  experimen- 
tal runs  and  was  defined  as: 

n     pud  _  md 
Re  =  - 


UAC 


where  d  was  the  height  of  the  channel,  u  was  the  dynamic 
viscosity  of  the  air  at  the  flow  exit  temperature  and  A 
was  the  cross  sectional  area  of  the  channel. 

For  all  experiments  in  the  curved  test  section,  the  Dean 
number  was  calculated  as : 


r^ 


De  =  Re 

Ri 


In  this  equation  d  was  again  the  height  of  the  channel  and 
Ri  was  the  radius  of  curvature  of  the  inner  convex  wall. 

To  verify  that  Q  .  ,  which  was  indirectly  measured  by  the 
temperature  rise  across  the  channel,  was  actually  the  con- 
vection heat  transfer  from  the  heated  Temsheet  surface,  an 
estimate  of  the  radiation  heat  transfer  was  made  for  each  run 
The  average  temperature  of  the  lower  plate  was  measured  by 
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eight  thermocouples  and  the  radiated  heat  transfer  (Q\)    was 

calculated  with  the  Stefan-Boltzman  law  of  radiation: 

u    u 
cr(T   -T  .  ) 
n..       wo    Wl 


RR 
where  a  is  the  Stefan-Boltzman  constant  and  R_.  represents 

K 

the  total  resistance  to  radiation  heat  transfer  between  the 
plates  as  given  by: 

RR  =  AT"  C  T~   +  ITT  -   13 

PL     wo     wi 

In  this  form  e    represented  the  emissivity  of  the  heated 
wo    r  J 

Temsheet  surface  and  z    •    was  the  emissivity  of  the  adiabatic 

wi  J 

Plexiglas   surface . 

Appendix  B  contains  a  representative  heat  balance  and 
sample  calculations  for  one  of  the  experimental  runs . 

B.   RESULTS 

Utilizing  the  expressions  described  earlier,  the  experi- 
mental data  for  buoyant  forces  aligned  with  the  main  direction 
of  the  fluid  flow  was  reduced  in  accordance  with  the  sample 
calculation  procedure  shown  in  Appendix  B.   Results  for  the 
major  variables  involved  are  provided  in  tabular  and  graphical 
form  on  the  following  pages .   Table  I  and  II  show  computed 
values  for  the  straight  and  the  curved  test  sections  respec- 
tively.  Figure  8  contains  a  plot  comparing  the  average 
Nusselt  number  versus  the  Reynolds  number  for  both  test 
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sections .   Representative  uncertainty  bands  have  been  indi- 
cated for  all  plotted  experimental  data  and  a  complete  error 
analysis  is  provided  in  Appendix  A. 
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TABLE  I 
SUMMARY  OF  RESULTS  FOR  STRAIGHT  TEST  SECTION 


Re 

h(W/m2°C) 

Nu 

989 

22.51 

5.321 

876 

22.21 

5  .298 

826 

21.79 

5  .198 

751 

21.75 

5.197 

689 

20  .72 

4.941 

632 

20  .95 

4.897 

560 

19  .20 

4.559 

488 

18  .79 

4.482 

438 

17.75 

4.234 

363 

16.99 

4.053 

313 

16.24 

3.874 

263 

14.30 

3  .512 
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TABLE  II 
SUMMARY  OF  RESULTS  FOR  CURVED  TEST  SECTION 


Re 

Dc 

h(W/m2°C) 

Nu 

990 

144 

27.81 

6  .694 

865 

126 

27.74 

6  .617 

803 

117 

25.06 

5  .977 

741 

108 

23.42 

5.588 

678 

99 

2  2.57 

5  .384 

654 

95 

22.43 

5  .351 

616 

90 

22.07 

5  .265 

579 

84 

21.67 

5  .169 

554 

81 

21.34 

5.091 

523 

76 

20.58 

4.875 

492 

72 

20.28 

4.838 

429 

63 

20  .20 

4.819 

404 

59 

20.31 

4  .844 

367 

54 

18.95 

4.515 

342 

50 

17  .76 

4.236 

306 

45 

15  .76 

3.759 

274 

40 

14.50 

3.460 

244 

36 

13.12 

3.130 
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In  all  cases  the  results  indicated  an  increase  in  the 
heat  transfer  rates  as  the  flow  velocity,  i.e.,  Reynolds 
number,  was  increased.   However,  when  comparing  the  average 
Nusselt  numbers  of  the  straight  and  curved  sections  two 
distinct  regions  became  evident.   At  the  highest  values  of 
Reynolds  number,  which  corresponded  to  the  highest  volu- 
metric flow  rates,  it  was  observed  that  the  values  of  the 
average  Nusselt  number  for  the  curved  test  section  were 
greater  than  those  in  the  straight  test  section.   At  the  low 
and  intermediate  values  of  Reynolds  number,  the  data  was 
less  conclusive  with  regards  to  the  comparative  values  of 
the  average  Nusselt  number  in  each  test  section.   This  region, 
between  Reynolds  numbers  of  200  and  400,  corresponded  to 
values  of  the  Dean  number  between  30  and  60.   As  Dean  discussed 
in  his  work  concerning  flow  instability  in  a  curved  channel 
[12],  a  small  disturbance  would  not  be  amplified  below  the 
critical  Dean  value  of  36.   Above  this  critical  value,  his 
linearized  prediction  indicated  the  disturbance  would  increas- 
ingly amplify  as  flow  velocity  increased. 

Consequently,  it  was  concluded  that  the  Taylor-Gortler 
vortices  are  only  initially  developing  in  the  low  to  inter- 
mediate Reynolds  number  region  and  have  not  yet  reached  suf- 
ficient proportion  to  have  any  significant  effect  on  the 
heat  transfer  rates  in  the  curved  section. 
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The  actual  presence  and  onset  of  the  Taylor-Gortler 
vortices  was  independently  verified  through  the  observation 
of  the  liquid  crystal  isotherm  distribution  on  the  Temsheet 
surface  in  the  curved  test  section.   At  the  smallest  values 
of  Reynolds  number,  the  liquid  crystals  were  in  essentially 
straight  bands  that  were  perpendicular  to  the  flow  direction 
(Figure  9).   A  similar  pattern  was  observed  in  the  straight 
test  section  for  the  entire  experimental  range  of  Reynolds 
numbers.   This  pattern  was  thought  to  indicate  that  the  flow 
was  in  essence  two-dimensional.   As  the  value  of  the  Reynolds 
number  was  increased,  liquid  crystal  bands  orientated  in  the 
streamwise  directions  became  visible  at  the  trailing  edge  of 
the  curved  test  section.   At  the  onset  of  these  streamwise 
bands  the  liquid  crystal  bands  at  the  leading  edge  remained 
orientated  perpendicular  to  the  flow  direction.   As  the 
Reynolds  number  was  further  increased,  the  bands  which  were 
orientated  in  the  streamwise  direction  eventually  became 
evident  over  the  entire  surface  of  the  curved  test  section 
(Figure  10).   These  streamwise  bands  would  indicate  a  transi- 
tion from  the  initial,  stable,  two-dimensional  flow  to  a  less 
stable,  three-dimensional  flow  triggered  by  the  continuing 
development  of  the  Taylor-Gortler  vortices. 

The  simultaneous  occurrence  of  the  streamwise  bands, 
explained  by  the  presence  of  the  longitudinal  Taylor-Gortler 
vortices,  and  the  increase  in  the  average  Nusselt  number 


46 


Figure  9 .  Distribution  of  Liquid  Crystal 

Color  Bands  in  Curved  Test  Section 
for  Re=250   (De=35) 
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1  1        / 

9     ^g^r  - ;                   ^  /--T 

.Figure  10.  Distribution  of  Liquid  Crystal 

Color  Bands  in  Curved  Test  Section 
for  Re=800   (De=117) 
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justified  the  conclusion  that  the  vortices  did  in  fact  en- 
hance heat  transfer  from  the  curved  test  section. 

Comparison  of  the  average  Nusselt  numbers  for  the 
straight  and  curved  test  section  indicated  that  values  for 

the  ratio  of  Nu   to  Nu   could  be  separated  into  two  cate- 

c      s 

gories.  For  Reynolds  numbers  between  400  and  800,  the  heat 
transfer  rate  was  increased  by  an  average  of  eleven  percent. 
For  Reynolds  numbers  above  800,  the  heat  transfer  rate  was 
increased  by  more  than  twenty-five  percent.   In  an  earlier 
investigation  of  the  boundary  layer  along  a  concave  wall, 
Kreith  [*+],  reported  that  heat  transfer  rates  along  concave 
walls  increased  from  twenty-five  percent  to  sixty  percent 

for  Reynolds  numbers,  based  on  the  hydraulic  diameter, 

M-        6 
between  10   and  10  .   This  would  suggest  that  the  two  data 

points  for  Reynolds  numbers  greater  than  800  in  the  curved 

test  section  could  represent  flow  that  has  tripped  and 

become  turbulent,  or  at  least  entered  a  laminar-turbulent 

transition.   One  possible  explanation  could  be  that  the 

Taylor-Gortler  vortices  have  started  to  break  down  at  this 

point  and  have  therefore  initiated  transition  to  the  turbulent 

flow  regime. 

Experimental  data  obtained  for  the  curved  test  section 

with  the  buoyant  forces  opposed  to  the  main  direction  of  the 

flow  indicated  a  minimal  difference  in  heat  transfer  rates 

from  the  curved  channel.   This  small  decrease  was  only 
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evident  at  very  low  volumetric  flow  rates  and  became  negli- 
gible for  Reynolds  numbers  greater  than  500. 
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V.   DISCUSSIONS  AND  CONCLUSIONS 

Investigation  of  the  measured  wall  and  air  flow  tempera- 
tures indicated  there  was  negligible  difference  between  the 
unheated  wall  temperature  and  the  fluid  bulk  temperature. 
In  addition,  energy  balances  such  as  the  one  detailed  in 
Appendix  B  provided  evidence  that  the  radiation  heat  transfer 
from  the  unheated  wall  was  minimal.   Based  on  this  information 
it  was  assumed  that  the  heat  convected  to  the  air  flow  was 
solely  from  the  heated  wall. 

As  previously  discussed,  the  high  aspect  ratio  of  the 
channel  provided  the  basis  from  which  the  initial  assumption 
was  made  to  model  the  experimental  apparatus  configuration  as 
infinite  parallel  plates .   The  experimental  data  was  found  to 
substantiate  this  assumption  as  the  computed  values  of  the 
average  Nusselt  number  in  the  straight  section  asymptotically 
approached  the  theoretical  value  of  5.385  for  parallel  plates 
with  one  wall  at  a  constant  heat  flux  and  the  opposite  wall 
adiabatic.  [28] 

With  the  formulation  of  the  problem  complete,  the  experi- 
mental data  could  now  be  compared  with  analytical  solutions 
of  the  same  problem;  i.e.,  laminar  flow  between  infinite 
parallel  plates  with  one  wall  at  constant  heat  flux  and  the 
opposite  wall  adiabatic.   In  particular,  the  experimental 
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results  were  compared  with  two  analytic  studies;  one  con- 
ducted by  McCuen,  et  al. [26],  for  infinite  parallel  plates, 
and  the  other  conducted  by  Wors</>e-Schmidt  [27]  for  tubes 
and  annular  passages.   In  the  case  of  the  annular  passages 
the  results  were  compared  to  the  numerical  predictions  for 
an  annulus  with  an  infinite  inner  and  outer  radius  of  curva- 
ture, thereby  approximating  parallel  plates. 

Before  comparing  the  experimental  data  with  the  studies 
mentioned  above,  it  would  be  advantageous  to  address  the 
initial  assumptions  and  review  the  nomenclature  common  to 
both  those  studies.   First,  the  following  restrictions  are 
applied  in  the  mathematical  formulation  of  the  problem: 

1.  the  velocity  profile  is  fully  established 

2.  the  fluid  transport  properties  and  density  are 
assumed  constant 

3 .  axial  heat  conduction  is  negligible 

4.  viscous  energy  dissipation  is  negligible 

5.  conditions  are  invariant  with  time  (steady  state) 
All  of  these  restrictions  have  been  accounted  for  in  the 
design  of  the  experimental  configuration  and/or  verified  in 
the  experimental  process.   Next,  a  dimensionless  axial  length 
coordinate  was  defined  as  follows: 


x 

X"    Oh  Pr  Re 
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where  x  was  the  length  of  the  straight  test  section  (29.2  cm), 
Dh  was  the  calculated  hydraulic  diameter  (1.239  cm),  and  Pr 
was  the  Prandtl  number  (0.70).   In  these  studies,  the 
Reynolds  number  (Re)  was  calculated  on  the  basis  of  the 
hydraulic  diameter.   The  values  of  x*  indicated  that  the 
development  of  the  experiments  would  most  likely  occur  in 
the  thermal  entrance  region.   Analysis  of  the  temperature 
distribution  on  the  heated  wall  versus  the  axial  length 
verified  the  flow  was  indeed  thermally  developing. 

In  comparing  the  experimental  results  to  those  reported 
by  McCuen,  et  al.,  it  was  necessary  to  integrate  their 
expression  of  the  local  Nusselt  number  for  small  values  of 
x*; 

Nu  = 


1/7  Q  • 

0.67095  x*1   -2(1  +  q^)x* 

^wo 


in  order  to  obtain  the  following  expression  for  the  average 
Nusselt  number. 


Nu  =  Q- ln  [1  "  0.67095(1  +  —  )X"    ] 

4(1+Iwi)x* 

^wo 

In  this  expression  Q  .  would  represent  the  heat  convected  to 

wi 

the  air  from  the  inner  Plexiglas  wall  and  Q    would  be  the 

&  xwo 

heat  convected  from  the  outer  Temsheet  wall.   The  average 
Nusselt  number  was  then  calculated  in  accordance  with  the 
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experimental  model  where  Q  .  was  equal  to  zero  and  Q  „  was 

r  Wl  wo 

synonymous  with  Q  •  • 

air 

The  experimental  data  obtained  for  the  straight  section 
was  then  recalculated  as  a  function  of  the  hydraulic  diameter 
and  both  the  experimental  and  analytical  results  were  plotted 
in  Figure  11.   Note  that  the  analytic  expression  for  the 
average  Nusselt  number  was  only  plotted  for  values  of  the 
Reynolds  number  that  corresponded  to  x*  less  than  0.0M-2  as 
this  represented  the  maximum  value  for  which  the  analytic 
solutions  would  be  valid.   Neglecting  those  experimental  data 
points  that  plotted  outside  this  valid  range,  there  was  good 
agreement  between  the  experimental  and  the  analytic  results. 

A  similar  process  was  followed  to  compare  the  experimental 
results  with  those  reported  by  Wars#e-Schmidt .   Again  the 
analytic  values  for  the  local  Nusselt  number  were  integrated 
to  obtain  analytic  values  for  the  average  Nusselt  number. 
In  this  case  the  experimental  and  analytic  values  of  average 
Nusselt  number  were  plotted  versus  the  non-dimensional  axial 
length  coordinate  x*  .   The  results  are  shown  in  Figure  12. 
Similarly,  there  was  good  agreement  between  the  experimental 
and  the  analytic  results . 

It  is  noted  that  in  both  cases  where  the  analytic  and 
experimental  results  were  compared,  the  experimental  data 
plotted  above  the  analytic  predictions.   This  difference  can 
be  explained,  at  least  in  part,  by  the  limitations  that  are 
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inherent  in  any  experimental  configuration.   In  this  case, 
the  presence  of  side  wall  effects  and  the  inability  to 
totally  account  for  all  the  heat  transfer  processes  and/or 
losses  combine  to  produce  a  slightly  inflated  value  of  the 
heat  transfer  rate.   Consequently,  the  average  Nusselt  number 
is  also  slightly  inflated. 

Irrespective  of  the  slight  differences  locally,  the  over- 
all trends  were  conclusive.   The  heat  transfer  rates  were 
shown  to  increase  with  increasing  values  of  Reynold  number. 

The  ratio  of  N    to  N    indicated  that  there  was  a  measureable 
uc     us 

increase  in  the  heat  transfer  rate  from  the  curved  section 
once  the  Taylor-Gortler  vortices  had  developed  to  sufficient 
strength.   Although  much  less  conclusive  from  this  experi- 
mental data,  it  appeared  the  the  Taylor-Gortler  vortices 
were  additionally  important  in  the  transition  from  laminar 
to  turbulent  flow. 

The  experimental  data  from  the  curved  test  section  when 
the  buoyant  forces  were  opposed  to  the  main  flow  direction 
was  essentially  inconclusive.   While  the  trend  was  qualita- 
tively what  had  been  expected,  that  is,  the  average  Nusselt 
number  for  the  buoyant  forces  aligned  with  the  flow  was 
greater  than  that  for  the  buoyant  forces  opposed  to  the  flow; 
the  difference  was  insignificant  in  light  of  the  experimental 
uncertainties . 
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VI.   RECOMMENDATIONS 

While  the  amount  of  work  remaining  in  the  experimental 
study  of  heat  transfer  in  a  curved  channel  has  been  reduced 
for  laminar  flows,  a  large  amount  still  remains  to  be 
accomplished.   Additional  experiments  should  be  conducted 
to  refine  the  correlation  for  the  average  Nusselt  number  as 
a  function  of  the  Dean  number.   Further  work  should  certainly 
be  considered  in  the  region  where  the  laminar  flow  begins  its 
transition  to  turbulent  flow.   Detailed  study  in  this  area 
could  provide  important  data  to  determine  the  manner  of 
relationship  that  exists  between  Taylor-Gortler  vortices  and 
laminar-turbulent  transition.   With  this  work  it  will  be 
necessary  to  design  and  construct  a  larger  experimental 
configuration  capable  of  higher  volumetric  flow  rates  and 
increased  AT  through  the  test  sections. 

Finally  it  would  be  advantageous  to  do  more  experimental 
work  in  the  general  area  of  liquid  crystal  thermography,  and 
more  specifically,  the  photography  of  these  crystals  and  their 
distribution . 
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APPENDIX  A 
ERROR  ANALYSIS 

The  uncertainties  for  the  major  variables  in  the  experi- 
ments were  calculated  in  accordance  with  the  method  described 
by  S.  Kline  and  F.  iMcClintoch  [31].   The  estimates  of  the 
uncertainty  in  the  measured  quantities  were  made  conserva- 
tively so  that  there  was  considerable  confidence  in  the  cal- 
culated uncertainties.   As  an  example,  the  calculation  of  the 
uncertainty  for  the  Reynolds  number  is  given  below.   The 
Reynolds  number  was  defined  by  the  equation: 

Re  = 


VI  •   A 

air   c 


and  the  uncertainty  was  calculated  as: 


"   m        d        u  .        A 


dRe 

Re    ▼  xlrTy     "  d  '     v  u  .  '  *  A 

Hair        c 


The  uncertainty  in  the  mass  flow  rate  was  determined  as 
follows : 


dm 
~m~ 


y(^)2 ;  ^  ;  ^2  +  ^p 


The  uncertainty  in  the  pressure,  the  gas  constant,  the  tem- 
perature, the  reading  of  the  volumetric  flow  rate,  the 
height  of  the  channel,  the  dynamic  viscosity  and  the  cross 
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sectional  area  of  the  channel  were  obtained  from  estimates  as 
follows : 


$f   =  0.00048,  ^  =  0.00002,  ^  =  0.00056,  ^  =  0.02000, 


dri  du  .  d  A 

^  -    0.02003,  — —  =  0.00223,  and  — r-£  =  0.02054 

yair  c 


The  uncertainty  in  the  mass  flow  rate  was: 


Si  =   V  (0. 00 048 ) 2+(0.00002)2+(0. 00055) 2+(0. 02000) 2  =  0.020 01 
m     T 

and  the  uncertainty  in  the  Reynolds  number  could  then  be 
calculated  as: 


—^  =V(0-02001)2+(0. 02003) 2+(0.00223)2+(0. 02054) 2= 0.03505 

Re    » 


The  values  of  the  uncertainties  for  all  variables  are 

listed  below: 

Quantity  Uncertainty 

A  0.02054 
c 

ApL  0.00266 

C   •  0.00415 
pair 

De  0.04042 

h  0.02149 

K  .  0.00038 
air 

Nu  0.02938 
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Quantity  Uncertainty 

m  0.02001 

Q  0.02000 

Q  .  0.03548 

R  0.03505 
e 

TD  0.00482 
B 

T -N  0.00257 

TOUT  °-00950 

T  0.00583 

TOUT-TIN  °-0290° 

AT=TT7-TD  0.00736 

W  B 

y  0.00223 

p  0.00074 
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APPENDIX  B 
SAMPLE  CALCULATIONS 

Figure  13  provides  a  representation  of  the  control  volume 
used  in  the  energy  balance  on  the  straight  test  section.   The 
sample  calculations  that  follow  indicate  the  procedure  used 
to  obtain  the  heat  transfer  components  required  for  the 
energy  balance  and  eventually  the  values  of  the  average 
Nusselt  number  as  a  function  of  the  Reynolds  number  for  the 
air  flow.   A  similar  energy  balance  and  procedure  was  used 
for  the  curved  test  section. 

ENERGY  BALANCE 
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Figure  13.  Energy  Balance  in  the  Straight  Test  Section 
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SAMPLE  CALCULATIONS 


A .   DATA 

TROOM 

Q 

V 
PR 

V 
u 

RPR 

T 
IN 

TINS1 

TINS2 

Twol 

Two  2 

Two  3 

Two4 

Two  5 

Two  6 

Two  7 

Two  8 

Twil 

Twi2 

Twi3 

Twi4 

Twi5 

Twi6 

23.0°C 
145% 

2.583  V 
30.008  V 
2.024  n 

0.892  mV  =  23.07  °C 
1.597  mV  =  39.514  °C 
1.339    mV      =       3  3.45    °C 

1.853  mV 

1.843  mV 

2.010  mV 

2.014  mV 

2 .144  mV 

2.141  mV 

2. 274  mV 

2.239  mV 

1.079  mV 

1.08  3  mV 

1.130  mV 

1.131  mV 
1.183  mV 
1.172  mV 


Two  =  5  0.47 


Twi  =  28.62 
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Twi7 

= 

1.219  mV 

Twi8 

= 

1.216  mV 

T 

= 

1.38  3  mV   = 

34,09  °C 

T 
EXIT 

= 

0.934  mV   = 

23.  53  °C   = 

:   74.35  °F 

APL 

= 

0.0742  m2 

A 
c 

= 

0.0016  m2 

d 

= 

0.  00635  m 

K  . 

air 

= 

0.02662  w/m 

°C 

Cp  . 
*air 

= 

1.0057   KJ/K 

g  °c 

air 

= 

1.983  x  10"5 

Kg/m* sec 

KINS 

= 

4.18  x  10"2 

W/m  °C 

AXINS 

= 

0  .00635  m 

WO 

= 

0.70 

£   . 

Wl 

= 

0.40 

F 

WO-Wl 

= 

1.0 

a 

= 

5.669  x  10"8 

W/m2  °K 

B.   TEMPERATURE  CALCULATIONS 

1.  Temperature  Difference  in  Insulation  (AT,^) 

iTINS  =  TINS1-TINS2  =  39.54-33.45  =  6.09  °C 

2.  Bulk  Temperature  (TR) 

t -  TIN  *  TOUT  _    23.07  +  34.09  .  9P  ,Q  or 
TB ^ 2 28'58   C 

3.  Mean  Temperature  Difference  (AT) 

AT  =  T    -  T_  =  50.47  -28.58  =  21.89  °C 
wo     B 
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C.   POWER  CALCULATIONS 

1.  Power  Supplied  (Qp) 

n     VPR  VH     (2.583X30.008)  _  ,fl  qn  ,7 
P  =  ~~*PR (2.024) 38'30  ^ 

2.  Heat  Lost  Through  Outer  Plate  (Q   ) 

AT 

axINS      _  6.09       

^Lo  "  AXINS/KINS*APL  '   ( 0. 00635 )/( (4. 18xl0"2)( 0.0742)) 


=  2.97  W 


3.   Heat  Radiated  (Q.  ) 


a.   Radiation  Resistance  (R„) 

— K 

1-e         ,         1-e  . 

r   =  wo  +  1 + 


R  Ant  £  n-nT       r  •  r\T^T       £      • 

PL    WO  PL       WO-W1  PL       Wl 


i    [JL  ♦  _i_  .  i]  =  2-929 


AnT         £  £      •  A^T 

PL        wo  wi  PL 


b.      Heat   Radiated    (Qr) 

4  4 

a(T         -T     .     ) 
Qr    =    wo wi_  T    =    OK 

RR 


(5.669xl0"8)(0.0742)(323.63l4-301.7814) 

2.929 


=    3.843    W 
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4,   Heat  Convected  to  Air  (Qair>) 


a.   Corrected  Volumetric  Flow  Rate  (q) 


/  tfytt    T°R) 

q   =   Q/0   =   Q/    t^ PRT 

1  STANDARD v       ; 


(0.314) (0.45)  /   534.04 

60  /     529.59 


3 
0  .0024   m    /sec 


b.       Density    (p ) 

P  1.013xlQ5 

?  ~    R   TEXIT    "    (287)(23. 53+273.15) 


=    1.189    Kg/m3 


c.   Mass  Flow  Rate  (m) 


m   =  qP  =  (0.0024) (1.189)  =  0.0028  Kg/sec 


d.   Heat  Convected  to  Air  (Q  •_) 

■ .^ — —    aix 


Qa.r  =  m  Cpair  (T0UT~TIN) 

=  (0.0028)(1.0057)(34. 09-23. 07) 

=  0.03103  KJ/sec 
=  31.03  W 
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D.  AVERAGE  HEAT  TRANSFER  COEFFICIENT  (h) 

r-    Qair      31.03 

n  "  ApL  AT      (0.0742M21.89) 

=  19.11  W/M2°C 

E.  AVERAGE  NUSSELT  NUMBER  (Nu) 


«-     hd       (19. 11)(Q. 00635)   _   u  „Q 

Nu  =  i7 — i —   =   7"n — nocco  S =    ^.559 

Kair  (0.02662) 


F.   REYNOLDS  NUMBER  (Re) 


m  d        (0.0028)(0. 00635) 
Re   = 


yair  Ac      (i:983x!0"5)(0.0016) 


=  560.39 
=  560 
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